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1 INTRODUCTION 



The nearest of all the giant radio galaxies, Centaurus A 
(NGC 5128) provides a unique opportunity to observe the 
dynamics and morphology of an active galaxy in detail 
across the electromagnetic spec trum. For a recent review 
on this remarkable object see llsraell (|l998l ). In its cen- 
tral regions, NGC 5128 exhibits a well recognized and 
optically-dark band of absorption across its nucleus. Images 
from the Spitzer Space Telescope (SST) with the Infrared 
Array Camera (IRAC) and Multiple Imaging Photometer 
for Spitzer (MIPS) cameras in the mid-infrared reveal a 



ABSTRACT 

We report on the results of spectroscopic mapping observations carried out in the nu- 
clear region of Centaurus A (NGC5128) over the 5.2 - 15 and 20-36/im spectral regions 
using the Infrared Spectrograph on the Spitzer Space Telescope. We have detected and 
mapped S(0), S(2), S(3), and S(5) pure rotational transition lines of molecular hydro- 
gen and emissions in the fine-structure transitions of [Sill] , [SIII] , [Fell] , [Felll] , [Aril] , 
[SIV], [Nell], [NeV] and [OIV]. The 500 pc bipolar dust shell discovered by Quillen et 
al.(2006) is even more clearly seen in the 11.3^m dust emission feature than previous 
broad band imaging. The pure rotational lines of molecular hydrogen other than the 
S(0) line are detected above the dusty disk and associated with the oval dust shell. 
The molecular hydrogen transitions indicate the presence of warm gas at tempera- 
tures 250-720K. The column density of the warm molecular hydrogen in the shell 
is N{H2) ^ 10^"cm^^ and similar to that estimated from the continuum dust shell 
surface brightness. The ratio of the dust emission features at T.T/um and 11.3/xm and 
the ratio of the [NeII](12.8^m) and II.3/zm dust emission feature are lower in the 500 
pc dust shell than in the star forming disk. The clearer shell morphology at 11.3/im, 
warm molecular hydrogen emission in the shell, and variation in line ratios in the 
shell compared to those in the disk, confirm spectroscopically that this shell is a sep- 
arate coherent entity and is unlikely to be a chance superposition of dust filaments. 
The physical conditions in the shell are most similar to Galactic supernova remnants 
where blast waves encounter molecular clouds. The lines requiring the highest level of 
ionization, [NeV](24.318/im) and [OIV](25.890^m), are detected 20-25"north-east and 
south-west of the nucleus and at position angles near the radio jet axis. Fine struc- 
ture line ratios and limits from this region suggest that the medium is low density 
and illuminated by a hard radiation field at low ionization parameter. These higher 
S molecular hydrogen pure rotational transitions are also particularly bright in the 
same region as the [OIV] and [NeV] emission. This suggests that the gas associated 
with the dust shell has been excited near the jet axis and is part of an ionization cone. 

3' long parallelogram shape (|Quillen et al.1 l2006bl ) that has 
been modeled as a scries of folds in a du sty warp ed thin 
disk (e.g.. Bland 1986; Bland c t ahl Il987l : .Nicholso n et all 
ll992l : ISparkelll996l : lQuillen et al.ll2006lj) . The dusty disk is 
also the site of ongoing star formation at a rate of about 
1 M(T) yr~^, (based on the infrared luminosity estimated by 
fToorJ"! Centaurus A hosts an active nucleus 



lEckart et aL"l99C 
(e.K., lWhvsong fc Anton ucci 2004.: iMirabel et al.a 999') that 
has b een recently studied using infrared spect ra from the 
SST. l|Weedman et al.ll2005l : [Goriian et al.]|2007l ). Its nucleus 
exhibits a strong silicate absorption feature and emission 
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from [NeV]. Here we do not study the nucleus but focus on 
structure exterior to it. 

The IRAC and MIPS images of Centaurus A have re- 
vealed another surprise, a bipola r shell-like structure 500 pc 
north and south of the nucleus (jOuillen et al.ll2006al ). This 
shell, seen for the first time with Spitzer imaging, could 
be the first extragalactic nuclear shell discovered in the in- 
frared. 

In this paper we present mapping of the central 500- 
600pc of Centaurus A done with the Infrared Spectrograph 
(IRS) on board the SST. The infrared spectral maps were 
obtained to test the possibility that the apparent dust shell 
is a separate coherent structure and not a chance superposi- 
tion of dust filaments. They were also obtained to search for 
possible interactions between the AGN and the interstellar 
medium either by illumination from the AGN or caused by 
jets or outflows. 

Based on the discussion bv llsrael (|l998l ). we adopt a 
distance to Cen A of 3.4 Mpc. At this distance, 1 kpc 
corresponds to I'on the sky. All positions reported in this 
manuscript are given with respect to Epoch 2000. 



2 IRS OBSERVATIONS 

Observations were obtained in spectral mapping mode in the 
Short-Low (SL;5-14/im) resolution and Long-High (LH; 19- 
38^m) resolution modules of the Infrared Spectrograph on 
board the SST. The spectral resolution of SL is i? ~ 50—100 
whereas that of LH is ii ~ 600. In spectral mapping mode 
the Spitzer space craft points on a raster of discrete steps, 
settling at each position before the integrations begin. Each 
map was accompanied by a 4 position offset observation 
taken with the same slit pattern and exposure sequence. 
The offset (rogue) positions were used for sky subtraction. 
The rogue sky position is at RA=13:25:15.7 (Epoch 2000) 
DEC=-42:55:13, approximately 6'north of the nucleus and 
sufficiently far away that little emission from the galaxy was 
seen in the previously observed broad band IRAC and MIPS 
images. 

2.1 SL observations 

We mapped in the short low (SL) module of the Infrared 
Spectrograph using slit overlaps of 25" x 1.85"using a 60s 
ramp. One data set covers 5 x 28 map positions (through 
the guaranteed time observer program; GTO), and the other 
two cover 7 x 30 map positions (through the guest observer 
program; GO). Both maps were observed by the SST on Mar 
13, 2007. Spectral cube assembly, including background sub- 
traction, bad pixel removal, spatial regridding and calibra- 
tion was done using the software package CUbe Builder for 
IRS Spectra Maps, (CUBISM^ (Smith et al.. ,2007b l pro- 
ducing a single spectral cube. 

2.2 LH observations 

We mapped in the long high spectral resolution (LH) module 
of the Infrared Spectrograph using slit overlaps of 10" x 4.5", 

^ http:/ /ssc. spitzer. caltech.edu/archanaly/contributed/cubism/ 



so that central positions in the map were covered by four slit 
positions. We used a 60s ramp, and each observation con- 
tained 2 cycles. For central map positions, this corresponds 
to a 480s on source integration time. Two data sets (guest 
observer; GO) were taken covering 6x20 map positions. The 
third (guaranteed time observer; GTO) data set covers 5x16 
map positions. All three maps were observed by the SST on 
Aug 1, 2006. 

Data reduction of the three LH spectral data sets in- 
volved bad pixel removal, spectral extraction, spatial regrid- 
ding, and application of the slit loss correction. For descrip- 
tions of the s o ftwar e procedures used see Appendix A by 
iNeufeld et al.l l|2006l ). The resulting three data cubes were 
then resampled onto the same spatial grid and averaged into 
a single cube, resulting in 1440s on source exposure times 
for central positions in the map. 

2.3 Continuum morphology in the spectral cubes 

We have extracted continuum images from both spectral 
cubes to check that they are consistent with previous broad 
band imaging an d to bring attention to features noted in 
previous studies lOuillen et al.1 l2006bl lah and that we will 
discuss below. 

In Figure [T] we show the morphology of the contirmum 
emission in 0.4/im wide bands centered at 25/im and 32/im 
from the long-high spectral cube. We chose these regions 
so as to be free from artifacts caused by order separation 
and emission lines. The morphology at 25/^m is consistent 
with the morphology seen in th e broad band MIPS 24//m 
image previously presented by (|Quillen et al.l [2006bf ) . The 
broad band of emission approximately 30" wide, and stretch- 
ing approximately along the east-west direction across the 
gala xy nucleus we refer to as 'the parallelogram' (see figure 
3 bv lOuillen et al.l [2006b). The parallelogram shape is due 
to folds in a thin dust y warped disk. The w arped disk is also 
rich in molecular gas (|Ouillen et al.|[l992t ) and so is the site 
of most active star formation present in the galaxy. 

In the SL spectral cube, broad and prominent dust emis- 
sion features are present. These dust emission features are 
generally attributed to polycyclic aromatic hydrocarbons 
(PAHs). These dominate the broad band e mission in the 
IRAC broad band images previously studied (|Quillen et al.l 
l 2006a '!^ except possibly at the nucleus. In Figure [2] we show 
the morphology of the peak of the 11.3/im dust emission fea- 
ture extracted from the SL spectral cube. This map also 
shows the parallelogram shape or star forming disk seen 
in the IRAC band images. Above and below the parallel- 
ogr am we can se e the sh ell-like bipolar structure described 
by lOuiUen et al.l (|2006al ) that was most prominent in the 
IRAC 8/im image. This shell has a major axis of ~ I'.l, a 
minor axi s of ~ ' .7 and position angle of approximately 
10° ( Qui llen et aLl lioOeaV The IRAC 8.0^m band image is 
dominated by the PAH dust emission feature at 7.7/im and 
so it is not surprising that the 11.3/im image shows simi- 
lar morphology. We show an image at 11.3/jm here because 
it most clearly reveals the dust shell. The morphology at 
7.7/im and other PAH emission features and their ratios will 
be discussed in more detail below. In the next section we will 
compare the morphology in line emission maps to that seen 
in continuum (Figure [1} and in the dust emission features. 

In the 25.0^m continuum image (Figure[T^) three peaks 
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in the emission are seen, 2 to the west of the nucleus, one 
to south-east. A weaker but fourth is located directly east 
of the nucleus. These are labeled with arro ws on the 24/im 
image in Figure 2 bv lQuillen et al. I l|2006ai ). We bring these 
to the attention of the reader as we extract spectra from 
these higher surface brightness regions to represent condi- 
tions typical of the star forming disk. We will show that 
spectra extracted from regions above or below the paral- 
lelogram region differ from those extracted from the star 
forming disk or parallelogram shaped region. 



3 EMISSION LINE MAPS 

Line emission images were created by summing emission 
near the line peak and subtracting nearby continuum. These 
images are shown in Figures [3] - [S] for lines present in the 
long wavelength spectral data cube and in Figures [8] - [11] 
for lines or dust emission features present in the SL spectral 
data cube. Line identifications and rest frame wavelengths 
are summarized in Table [1] 

3.1 Morphology associated with the star forming 
warped disk 

The lines [Sill] (34.8^im) and [SIII](33.5^im) (as displayed in 
Figure are strong in the parallelogram or star forming 
disk, as would be expected from the presence of HII regions 
in this disk. Peaks in these maps are approximately coin- 
cident with peaks seen in the nearby continuum at 32/im 
(Figure [Da). The [FeII](26.0/im) and [FeIII](22.93^m) lines 
are also detected in the same region with [Felll] weaker than 
[Fell] (Figure gj. Molecular hydrogen in the H2S(0)J = 2- 
0(28. 2/im) line is also detected associated with star forma- 
tion in the warped disk or parallelogram (Figure [5} though 
the morphology in this line differs from that seen in the 
[Sill], [SIII], and [Fell] lines. The northwestern and south- 
eastern ridges of the parallelogram are stronger than the 
southwestern and northern ridges. The northwestern and 
southeastern ridges also seem to be brighter in the [Fell] 
map. 

The emission lines present in the long wavelength spec- 
tral cube are clearly associated with the disk, however there 
are differences between these maps and the nearby contin- 
uum. There is more structure seen in the [SIII] map than 
present in nearby continuum at 32/xm; peaks in the disk are 
brighter than the surrounding omission in the parallelogram 
(compare Figure [3] to [TJd). Quillon ct al. (2006a) suggested 
that the 4 peaks in the broad band 24/im images were due 
to a superposition of the shell and the star forming disk. 
However the emission lines at long wavelengths are not de- 
tected in the vicinity of the shell-like structure but do display 
the four emission peaks (particularly the [SIII] line image). 
Therefore that explanation for the 4 peaks seen at 24/im is 
probably incorrect. Those peaks were not predicted b y the 
warped disk model presented by iQuillen et al.l l2006bl . The 
peaks are symmetrical about the nucleus suggesting that 
they are due to projection of a warped disk and that im- 
provements in the warp model might succeed in accounting 
for them. If this were the case then the size of the peaks 
should be similar in the continuum and line emission. Un- 
fortunately the peaks are more prominent in some of the 



lines (e.g., [SIII]) than the nearby continuum. Spiral struc- 
ture in the warped disk could cause particular regions in 
the disk to exhibit higher levels of star formation and could 
maintain the reflective symmetry. If this were the case then 
variations in the observed morphology between lines and 
continuum would require a model in which the strength of 
spiral structure depends on wavelength. 

3.2 [Nell] emission associated with the dust-shell 

The [SIII], [Sill] and [Fell] emission are primarily confined 
to the warped disk (or parallelogram shape). However the 
[NeII](12.8/im) image (see Figure [TT]) does show faint levels 
of emission (at the level of ~ 10~^erg cm~^s~^SR~^) associ- 
ated with the dust shell. This line is blended with the nearby 
12.7/im dust emission feature, nevertheless the shape of the 
spectrum at different locations suggests that [Nell] dom- 
inates the dust emission and is present above and below 
the galaxy disk. If so then the [SIII], [Sill] and [Fell] emis- 
sion associated with the dust shell is just too faint to see in 
the long wavelength spectral cube. Deeper observations at 
longer wavelengths or higher spectral resolution short wave- 
length infrared observations would be needed to confirm the 
presence of ionized gas in the shell. A spectroscopic study 
at visible band wavelengths would be quite sensitive as ex- 
tinction is probably not high in the shell region above and 
below the parallelogram. 



3.3 Morphology associated with higher ionization 
lines, [OIV] and [NeV] 

Lines requiring higher ionization, [OIV](25.9pim) and 
[NeV](24.3/im), are also detected in the LH spectral cube 
and they exhibit different morphology than the lower ion- 
ization lines. Emission from [OIV] and [NeV] is strong along 
position angles of ^ 40° and —120° from the nucleus and 
is resolved at locations outside the nucleus; see Figure [6] 
These position angles are within 15° of the jet axis, though 
the jet is seen prominently only on the north-east side in ra- 
dio emission. The radi o and X-ray jets in this region have a 
position angle of 55° ijBurns et al.lll983l : iKraft et aLlbOOOD . 
In figure [7] we show contours of radio continuum at 5 GHz 
bv lHardcastle et al.l l|2006l ) overlayed on the [OIV](25.9/im) 
emission map (shown in grayscale). We find that the position 
angle of the [OIV] differs from the jet axis by 5 — 15°. The 
different between the axis of the [OIV] emission and radio 
jet axis is perhaps not surprising. Observations of Seyfert 
galaxies have found that the jet axis can differ from that of 
an ionization cone (NGC 1068, for example). 

High angular resolution observations of the central few 
arcseconds of Centaurus A have revealed elongated or cone- 
shaped emission in Paa, Pa/3 and [FeII](1.26/im) north of the 
nucleus that has been interpreted in terms of a rotating disk, 
rather than in terms of an ionization cone (jSchreier et al.l 
ll998l : lKrai novic et al.ll2005l ). The jet axis lies on the west- 
ern side of the possible ionization cone seen in [FeII](1.26). 
The [OIV] and [NeV] extension to the north of the galaxy 
nucleus are approximately oriented in the s ame direction as 
the co ne seen in [FeII](1.26/im) shown bv iKrainovic et "al] 
(I 2OO5I ). As there is no cone feature seen in the near-infrared 
[Foil] observations south of the nucleus (likely because the 
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band of absorption caused by the warped disk) it is more 
difficult to compare the cone to the [OIV] emission south- 
west of the nucleus. The [OIV] emission north and south of 
the nucleus are not 180° apart and the emission south-west 
of the nucleus would not overlay the [Fell] (1.26) possibly 
ionization cone if rotated by 180°. 

It is interesting to reinspect the morphology of the lower 
ionization line maps at the position of the peak in [OIV] 
south west of the nucleus. An extension in the contours 
of the [SIII], [Sill], [Fell] [Felll] and H2S(0) emission line 
brightness is seen at the same location as the peak south- 
west of the nucleus in the [OIV] map. The ratio of the surface 
brightness in the extension compared to that in the disk is 
higher in the [Felll] line compared to the [Fell] line. This is 
not surprising as we know these extensions contain higher 
ionization species such as [NeV]. 

The spectral resolution in the LH spectral cube (a spec- 
tral resolution of 7? = ^ ~ 600 corresponds to a velocity 
resolution of Av = 500km/s) is sufficiently high that we can 
just barely see the galactic disk's rotation in the brighter 
lines (e.g., [SIII], [Sill] and H2S(0)), with the western and 
eastern sides of the disk redshifted or blueshifted, respec- 
tively, compared to the system velocity. We inspected the 
channels near the line centers and searched for evidence of 
structure in the radial velocity distribution. We found no 
clear evidence for a significant (greater than ±200km/s) red 
or blue shift in the [OIV] emission extensions compared to 
the system velocity. We also saw no clear evidence for an 
increase in line width. The line width (FWHM) must be 
smaller than the galactic rotation full width or < 400km/s. 
Near the [OIV] peak south-west of the nucleus the extensions 
in the contours of the [SIII], [Sill], and H2S(0) also appear to 
have velocities within 200km/s of the system velocity. Un- 
fortunately the spectral resolution is just not high enough 
to be certain about these line widths or allow a search for 
structure in the radial velocity distribution. 

3.4 Emission from molecular hydrogen 

In the SL spectral cube we detect four pure rotational molec- 
ular hydrogen lines, the S(7)J = 7-5, the S(5)J = 5-3, the 
S(3) J = 3-1 and the S(2) J = 2-0 lines. We show the contin- 
uum subtracted line emission maps for the H2S(5)(6.909/im) 
and H2S(3)(9.665^m) lines in FigureH The S(7) is weak and 
the S(2) strongly blended with dust emission features so we 
have only presented line emission maps from the other two 
lines. The S(4) line at 8.0251/im was not detected, probably 
due to the bright dust emission feature at 7.7/im at its ex- 
pected wavelength and because the even number transitions 
are singlets rather than triplets and so have one third the 
number of states. The S(l) line at 17.03/im was not detected 
because it lies outside the spectral range covered by the SL 
and LH modules. To observe this line we would have needed 
observations with the SH module of the Infrared Spectro- 
graph. 

The morphology of the higher S(3) and S(5) rotational 
molecular hydrogen lines differs from that we see at longer 
wavelengths in the S(0) line at 28/xm (compare Figure [5] to 
Figure [9|. The S(0) line is primarily associated with the 
star forming disk, however the S(3) and S(5) lines are seen 
above the disk. An overlay of the H2S(3)(9.665/xm) line and 
[OIV](25.89/im) in Figure [9^ shows that the peak in the 



[OIV] south-west of the nucleus coincides with the peak 
H2S(3) surface brightness in the same region. An overlay 
shown in Figure [SJ) of the S(3) line with the dust emission in 
the 11.3/i PAH emission feature (previously shown in Figure 
[2| shows that the warmer molecular hydrogen gas (traced 
by the higher S pure rotational transitions) is near the dust 
shell. The simplest explanation for the morphology is that 
warmer molecular hydrogen is associated with the dust shell 
but that there is additional excitation of the shell material 
near the jet axis. Constraints on the column depth and tem- 
perature of the molecular hydrogen are discussed below in 
section 5. 



3.5 Colors in the dust-shell compared to the 
warped disk 

In Figure [10] and Figure [11] we show maps of continuum 
emission at the peak of the 7.7, and 8.6/im dust emission 
features and the [NeII]12.81/im line. We did not subtract 
continuum from these maps as the [Nell] line and dust emis- 
sion feature are 3-8 times brighter than the continuum. As 
expected, the [Nell] line and dust emission features primar- 
ily trace the star forming regions in the warped disk. Fainter 
levels of emission arise from the dust shell above and below 
the disk and from the region that is bright in [OIV] emis- 
sion (~25"south west of the nucleus). In Figure [121 we show 
ratios of the 7.7^m and 11.3pim surface brightness and the 
ratio of the [NeII]12.8 and 11.3/im surface brightness. Both 
of these ratios exhibit variations across the galaxy, with the 
lowest values in the dust shell and higher values present in 
the star forming disk. Subtraction of the continuum made 
little difference in the morphology of the ratio maps because 
the lines and dust emission features are over twice as strong 

as t he continuum . 

[Smith et all (|2007[ ) found a correlation between the 
7.7/11.3/im ratio of dust emission complexes and the 
[NeIII](15.6)/[NeII](12.8/im) ratio in star forming galaxies. 
They also noted that the galaxies with low 7.7/11.3 com- 
plex ratios were more likely to be classified as LINERS 
or Seyferts. They suggested that the 7.7/im/11.3/jm ratio 
is an indicator of the hardness of the radiation field. The 
[NeIII](15.6/im) line lies is outside the wavelength coverage 
of the SL module so we cannot use that ratio to probe the ra- 
diation field above the disk. At th e nucleus in Gen A, thi s line 
was detected but was not strong ([Weedman et al.ll2005^ . We 
do have higher ionization tracers in the the [OIV](25.9/xm) 
and [NeV](24.3/im) lines. However, the morphology in those 
lines differs from that we see in the the 7.7 to 11.3/im ra- 
tio map. The [OIV] line is bright near the jet axis whereas 
whereas the 7.7/11.3 ratio is low in the dustshell (see Figure 
I12p . The lack of change in the 7.7/11.3 ratio in the region 
containing [OIV]25.9 and [NeV]24.3 emission near the jet 
axis makes it difficult to interpret variations in the 7.7/11.3 
ratio purely in terms of hardness of radiation field. Since 
the 7.7/11.3 ratio drops in the vicinity of the dust shell and 
this is where we detect warm molecular gas we might in- 
stead consider models in which PAH emission is affected by 
the presence of warm gas or the physical process that is re- 
sponsible for heating the warm gas. We could also consider 
models where this ratio is affected by the ambient radiation 
field strength. 
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4 SPECTRA 

We have extracted spectra at various positions in the star 
forming disk from the LH spectral cube and these are shown 
in figure [Tib . Spectra were extracted from 7"square regions 
centered at 5 regions. 4 of these correspond to peaks evi- 
dent in Figure [3] excluding the nucleus. From bottom to top 
the spectra shown in Figure 114b were extracted from re- 
gions centered at RA=13:25:25.8, DEC=-43 00 52 (peak to 
the west and north of the nucleus), RA=13:25:25.9, DEC=- 
43:01:09 (peak west of the nucleus) , RA=13:25:29.5, DEC=- 
43:01:27 (south-east of the nucleus), RA=13 25 29.5, DEC=- 
43:01:10 (peak east of the nucleus), and RA=13:25:28.6, 
DEC=-43:01:23 (peak south-east of the nucleus). 

We also extract an LH spectrum along the jet axis at 
the peak shown in Fig |14b . This spectrum was extracted 
over a ll"square region centered at RA=13:25:25.7, DEC=- 
43:01:24, at the peak in the [OIV] emission south-west of the 
nucleus (see figure [6} . The lines we have detected from these 
spectra are listed in Table [1] Line fiuxes measured from this 
region are listed in Tables [2] and [S] 

We see only minor differences in different regions of 
the parallelogram (see Figure [14)1 . However, along the jet 
axis there are deviations in the spectra, in the relative 
strengths of the lines requiring higher levels of ionization, 
[OIV](25.9/im) and [NeV](24.3^im), compared to those re- 
quiring lower levels (e.g., [SiH](34.8^m) ,[SIII](33.5/im)). To 
produce the emitting ion from the preceding ionization state, 
[OIV] requires an ionization energy of 54.0ev where as [NeV] 
requires 97. lev. A harder UV radiation field most naturally 
accounts for these species. 

In the SL spectral cube we extracted a spectrum from 
the same bright region along the jet axis south-west of the 
nucleus as that shown in Figure [T4b and from a representa- 
tive region in the parallelogram representing the star form- 
ing disk. The spectrum representing the star forming disk 
was extracted from a position centered at RA=13:25:25.9, 
DEC=-43:01:09. The location is the same as that for the 
second from bottom spectrum shown in Figure [141 The two 
spectra extracted from the SL spectral cube are shown in 
Figure 1131 

As can be seen from Figure flSb the star forming region 
exhibits little molecular hydrogen emission in the higher S 
rotational lines. However, the S(2), S(3), and S(5) lines are 
detected near the jet axis and above the disk near the dust 
shell (see Figure [T5bl . The S(7) line is barely detected at 
5.51lAtm. The S(4) line at 8.0241 is not detected likely be- 
cause it has been overpowered by the dust emission feature 
at 7.7/im. Likewise the S(6) line at 6.1086/im is not detected 
as it would have been overpowered by the dust emission fea- 
ture at 6.2/im. The S(5) line is detected above the disk but 
not in the star forming region at 6.909/im (Figure [T3b l . At 
almost the same wavelength the [Aril] line at 6.985/xm is 
seen in the star forming disk. The spectra are sufficiently 
high spectral resolution that we can differ between the two 
lines. A comparison between the two spectra also shows that 
the relative heights of the 7.7 and 11.3/im dust emission fea- 
tures differs. These two dust emission features are approxi- 
mately the same height in the star forming disk. However the 
11.3/im emission feature is twice as bright as the 7.7/im dust 
emission feature in the dust shell. The [NeII](12.81/im) and 
dust emission feature at 11.3^m are also similar in height 



in the star forming region but the 11.3/im is almost twice 
as bright as [Nell] above the disk. These changes are re- 
flected in the ratio maps shown in Figure [T2] and discussed 
in section 3. 

There is a weak feature near 14.2/xm in both spectra 
from the star forming disk and jet region (Figure I13b .b). 
The shape of the spectra are similar in shape at 14.2^m 
suggesting that this feature is a dust emission feature rather 
than due to emission from [NeV](14.3/im). A higher spec- 
tral resolution spectrum would be required to detect this 
line. Weak [NeV](14.3/j m) emission was detect ed in the SH 
spectrum of the nucleus l|Weedman et al.|[2005l ). 



5 LINE RATIOS 

We have measured line fluxes from the spectra shown in 
Figures 114b and 113b which were extracted at the location 
of the [OIV] peak south-west of the nucleus. Table [2] lists 
the fluxes of the emission lines from ions and Table O lists 
the fluxes of the pure rotational molecular hydrogen transi- 
tions. Line fluxes were measured from the 3 spectral cubes in 
the same aperture using continuum fitting and subtraction 
routines available in CUBISM. We have used CUBISM to 
measure line fluxes for the LH spectral cube even though the 
images (displayed in our flgures) show the spectral cube re- 
duced with the software described by Noufold ot al. (200f|). 
This spectral cube had fewer artifacts than the one produced 
by CUBISM, however the artifacts in the CUBISM reduced 
LH spectral cube should not have affected the measured line 
fluxes as we were able to choose where to measure the con- 
tinuum and line emission. Also the same software package 
was used to measure all line fluxes to check and ensure con- 
sistency in the calibration between the spectral cubes. 

5.1 Molecular hydrogen emission 

Pure rotational transitions of molecular hydrogen in the 
mid-infrared (between 3 and 20/im) are a major coolant 
of of warm gas in the temperature range 100-lOOOK 
(jNeufeld fc Kaufman|[l993 ). Higher S pure rotation molec- 
ular hydrogen lines have previously be en detected in ex- 
tragalactic objects. iHigdon et all l|2006l ) detected pure ro- 
tational lines from ULIRGS i ncluding the S(7) line in a few 
objects. jPanuzzo et al.l (|2007l ) detected these lines in the el- 
liptical galaxy NGC 4435 that is experiencing an off-center 
encounter with a nearby galaxy. This galaxy hosts a dusty 
disk as seen from visible band HST images and ongoing 
star fo rmation at modest rate of ~ O.IM© /vr. lRoussel et aTl 
^200i ) observed warm gas in the centers of nearby galaxies 
with LINERs and Seyferts having warmer H2 emitting in the 
excited states or higher S transitions. The above studies did 
not resolve the H2 emission or did not study its morphol- 
ogy. Higher angular resolution near infrared studies have 
found emission from molecular hydrogen in vibrational tran- 
sitions in active galaxies. In some Seyf erts this emission is 
resolved on l OOpc or larger scales (e.g., lOuillen et al. I ll999l : 
IPavies et al.ll2005l ) i n others rotation is see n near the mas- 
sive black hole fe.g.. lHicks fc MalkanllioOTi ). 

The flux of a pure rotational transition can be written 
F — ^ANjj^ where h is Planck's constant, c is the speed of 
light, A is the wavelength of the transition, A is the Einstein 
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A-coefEcient, Q is the solid angle and Nj is the column 
density of the initial J quantum state. Since we are consider 
pure rotational S transitions the rotational quantum number 
change is A J — —2, the vibrational quantum number change 
is Ai/ = 0, and v = Q. 

Assuming local thermal equilibrium, the fraction 
of molecules in the J— th state is N.j/N[H2) = 
g.jexp{-E,j/kT)/Z{T) where Z{T) is the partition func- 
tion. Here gj is the statistical weight of the state or 2 J + 1 
times 1 or 3 for the even and odd J states, respectively. To 
estimate the temperature of the warm molecular hydrogen 
we make an excitation diagram plotting the energy of the 
initial state Ej/k versus logio(A'^j/<7./)+ a constant. The ex- 
citation diagram is shown in Figure [151 On such a plot a sin- 
gle temperature gas with temperature T would have points 
with slope equal to log^o(e) /T. We have chosen to plot using 
logj^Q rather than the natural log so that the y-axis gives an 
estimate for the column density of the states. 

Figure [15] shows that a single temperature model does 
not fit all data points. This is not unexpected as other extra- 
galactic o bjects also require multi-temperature models (e.g.. 



Higdon et al.,2q06 : Panuzzo et al. 2007 ; R oussel et al. 2007. ; 
Ogle et al.ll2007l ). The points lie along a line for initial states 



higher than J = 2 implying that the higher J rotational level 
populations are close to LTE. The S(2) and higher lines are 
fit with a temperature of T2 = 720K. The S(0) and S(2) line 
are connected with a slope corresponding to a temperature 
of Ti ~ 260K. Had we subtracted the warmer component, 
a somewhat lower temperature would have been measured. 
Unfortunately our spectra do not cover the S(l) line. This 
would have helped in making a 2 component model. It is un- 
likely that there is a strong deviation from thermalization of 
ortho levels with para levels as the S(2) point lies along the 
same line as the S(3), S(5) and S(7) points in the excitation 
diagram. 

The column density corresponding to the above excita- 
tion model is N{H2) ^ lO^^cm"^ for the cooler molecular 
component (Ti — 260K) and ^ 10^^cm~^ for the warm 
component {T2 — 720K). The column densities are simi- 
lar to those measur ed in the nuclei of nearby galaxies by 
iRoussel et all (|2007l ) and in supern ova remnants where blast 
waves encounter molecular clouds (|Neufeld et al.|[2007l ') . The 
temperatures of the two compone nts we measure are hi gher 
than most of those measured by I Roussel et al.l (|2007l ) for 
most nearby galaxy centers excepting the LINERS and 
Seyfert galaxies, bu t lower than the supernova remnants 
l|Neufeld et al' '200?) and ULIRGs with two temperature 



models |Higdori et al. 2006). The fraction. 



lo- 



in the 



warmer component (that at 720K compared to that at 260K) 
is higher than that measured for the ULIRGS with 2 compo- 
nent tempe rature models and rn ost but not all the galaxies 
studied by iRoussel et all l|2007h but similar to that mea- 
sured from the supernova remnants (|Neufeld et al.l |2007| ) . 
The physical conditions in Cen A's dus tshell are most closel y 
matched by the supernovas studied bv lNeufeld et al.l(|2007^ . 

As we have discussed from the images, the molecular 
hydrogen emission in the higher S lines arises from the dust 
shell and along the jet axis, though that in S(0) is associated 
with the disk and so is probably from with photodissasocia- 
tion regions (PDRs) associated with star formation. We did 
not clearly detect the higher S lines in the disk suggesting 
that the molecular hydrogen at the hotter T2 = 720K tem- 



perature is primarily excited in the vicinity of the dust shell 
or along the jet axis. The column density of cooler molec- 
ular hyd rogen is similar to that estimated in the dustshell 
by Ouillen et al.l (|2006al ) from the shell surface brightness. 
This suggests that gas associated with the dust shell has 
been heated, in the shell and more so near the jet axis. 
There is probably more than one excitation process as [OIV] 
and [NeV] emission are detected primarily along the jet axis 
whereas the higher S pure rotational molecular hydrogen 
lines are detected there and in the vicinity of the dust shell. 



5.2 Excitation near the jet axis 

Nebular emission lines can constrain the the spectrum of the 
energy powering the emission. Our spectra allow us to make 
use of the following line ratios: [OIV]25.9/[SIII]33.5Atm, 
and [OIV]25.9/[NeII]12.8Atm. We do not detect [NeV] 
at 14.3/im so we will discuss limits on the ratios 
[NeV]14.3/[NeII]12.8^m and [NeV]14.3/[NeV]24.3^m. 

We measure [OIV]25.9/[SIII]33.5Mm ~ 1.0 at the [OIV] 
peak. This ratio is most similar to those of the AGNs con- 
sidered by (iGenzel et al.l Il998l) (see their figure 3) . How- 
ever we have failed to detect a significant [NeV] line at 
14.3/im so we estimate a limit on [NeV](14.3)/[NeII](12.8) 
< 0.05. Ratios this low are more similar to starbursts than 
AGNs. We note that the ratio we measure in this region 
[OIV]25.9/[NeII]12.8 ~ 0.17 is lower than expected from 
the [OIV]25.9/[SIII]33.5 ratio using Genzel et al.'s correc- 
tive factor of 1.7. Our estimated [OIV]25.9/[NeII]12.8 ~ 0.2 
ratio places the emission along the jet axis more similar to 
ULIRGs than AGNs. [OIV]2 5 9/ [Nell] 12.8 is exp ected to be 
approximately 1 for AGNs (Genzel et al.ll 19981 . Figure 3a 
and Figure 7 bv lSturm et al.(l2002D . Our low estimated ra- 
tios of [NeV]14.3/[NeII]12.8Mm and [OIV]25.9/[NeII]12.8Mm 
compared to those expected in AGNs may in part be due 
to the blending of the [Nell] 12.8 emission with dust emis- 
sion features which may have caused us to overestimate the 
[NeII](12.8/xm) line strength. 

We don't detect [NeV](14.3^m) or [NeVI] at 7.6^m and 
we do detect [NeV](24.3/im). By comparing upper limits of 
ratios of these li nes with [OIV] (25. 9) and [NeV] (24.3/im) 
and Figure 4 by ISturm et al.l l|2002l ) (based on models by 
ISpinoglio et al.l 2000l ) we infer the emitting region is likely 
to have low ionization parameter U < 10~^. Our estimated 
upper limit in the ratio of [NeV] 14.3/ 24.3 of < 3 in the jet re - 
gion and a comparison of Figure 3 bv lAlexander et al.1 ll999D 
suggests that the electron density is also low Ue < lO^cm"^. 
The detection of [NeV] suggests that a hard radiation field 
is illuminating the dust shell near the jet axis. However 
we have failed to detect [NeV](14.3/im) or [NeVI](7.65/im) 
suggesting that the illuminated region is low density and 
has a low ionization parameter. The difference in the phys- 
ical conditions compared to those expected in narrow line 
regions might account for the discrepancy between the 
[OIV]25.9/[NeII]12.8 and [OIV]25.9/[SIII]33.5 ratios com- 
pared to those seen in AGNs. 

We consider the possibility that the low electron den- 
sity and ionization parameter are consistent with illumina- 
tion of a region distant from the nucleus by hard radiation 
from the central AGN. The position of the [OIV](25.9/im) 
peak corresponds to a distance d ~ 400 pc from the nu- 
cleus. The active nucleus is estimated to have a bolomet- 
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ric lu minosity of Lboi ~ 10*^ erg/s (jWhysong fc Antonuccil 
|2004 ). The ionization parameter, U, for isotropic emission 
is [/ = 4^^„g where c is the speed of light, Q is the ioniz- 
ing photon emission rate (photons per second) and n is the 
gas density (in cm~'^). Assuming a mean energy for ioniz- 
ing photons of lOOev (the estimated energy of the UV bump 
and hard enough to produce NeV from the preceding ioniza- 
tion state) at the estimated bolometric luminosity emitted 
isotropically into Att radians. 



U ~ 0.001cm" 



I^10«erg s-i J \^400pi 



( " V 

Vl02cm-»y 



This estimate illustrates that a low ionization parameter, 
low density medium (such as we infer from the line ratios) 
is the only situation that we might expect at a distance of 
400 pc from the AGN. 

We now ask: can the AGN nucleus provide enough 
UV photons to account for the flux of the [OIV] emission 
25"from the nucleus? We summed the [OIV](25.9^m) flux 
in a 15"wide region centered on the [OIV] emission peak, es- 
timating a flux of 10"^'^ erg cm"'^ s~^ in the line. This flux 
corresponds to a luminosity in [OIV] of 10^* erg/s. We can 
compare th e [OIV] lumi n osity to the mid-infrared luminos- 
ity as have ISturm et ahl (|2002l ) for other objects. Since the 
[OIV] flux is illuminated only in a small solid angle from the 
nucleus we must first estimate the fraction of UV luminosity 
that could be causing the [OIV] emission. The [OIV] emis- 
sion comes from a region that is about 10"wide. As viewed 
from the position of the nucleus, this corresponds to a solid 
angle of (20°)^ or ~ 1% of 47r steradians. Thus we compare 
the [OIV] luminosity to 1% of the mid-infrared luminosity. 
The mid-infrared luminosity of the AGN is ~ 10*^ erg s~^ 
IWhvsong fc Antonuccil |2004| ). Consequently we estimate 
that the ratio of the [OIV] to 1% of the mid-Infrared lu- 
minosity is ~ 10~^. This rati o is similar to that s een for the 
brighter Seyferts studied by ISturm et al.l (|2002l ) (see their 
Figure 12). We conclude that as long as they are not ab- 
sorbed before they run into the dust shell, UV photons from 
the AGN can account for the [OIV] and [NeV] emission near 
the jet axis. An additional power source (like shocks asso- 
ciated with the jet) is not necessarily required to account 
for the high ionization species seen near the jet axis but 
distant from the nucleus. If UV photons from the AGN are 
responsible for the high ionization species then the column 
depth along the sight line between the AGN and dust shell 
must be low. If a blast wave is responsible fo r sweep ing up 
the dust shell, as proposed bv lOuillen et al.l l|2006al ). then 
it could also be responsible for evacuating material between 
the AGN and the dust sheU. 



6 SUMMARY AND CONCLUSION 

We have carried out a spectroscopic study of the central 2 
arcminutes of Centaurus A using short low and long high 
spectral modules of the Infrared Spectrograph on board the 
Spitzer Space Telescope. Most of the emission lines detected 
in the spectral cubes (e.g., [SIII](33.5^im), [Sill] (34.8), 
[FeII](26.0), [FeIII](23.9), [ArII](6.98) and H2S(0)(28.2)) 
and dust emission features primarily trace regions of star 
formation in the warped disk. 

Our previous study based on IRAC and MIPS imaging 



suggested that Centaurus A hosts an oval or bipolar dust- 
shell at a distance approximately 500 pc from the nucleus 
seen above and below the warped disk. This dust shell, if 
confirmed, would be the first extragalactic shell to be dis- 
covered in the infrared. Here we see the dust shell even more 
clearly and prominently in the 11.3/xm dust emission fea- 
ture than we saw previously in the broad band IRAC im- 
ages. We have found variations in the dust emission feature 
7.7/im/11.3/im ratio and dust 11.3/im/[NeII](12.8/i) ratio, 
with the oval dust shell having the lowest ratios compared 
to the star forming disk. The clearer shell morphology at 
11.3/im than previously seen in broad band images, the asso- 
ciation of the molecular hydrogen emission in the shell, and 
the variation in line ratios in the shell compared to those 
in the disk, con firm spectroscopically that the shell discov- 
ered previously (|Quillen et al.ll2006a ) is a separate coherent 
entity and is unlikely to be a chance superposition of dust 
filaments. 

We find evidence for higher ionization species line emis- 
sion in [NeV](24.3^m) and [OIV](25.9/im) near the jet axis. 
Emission in these two lines is seen both north-east and 
south-west of the nucleus along position angles ~ 40° and 
~ —120°. These angles are similar to but not exactly the 
same as the jet axis at 55° as seen at 5GHz. Outside the 
nucleus, the peak surface brightness in these lines is 25"or 
400 pc south-west of the nucleus. Emission line ratios and 
limits at the location of the [OIV] peak suggest that the 
emitting region is at low ionization parameter, U < lO"'^, 
and has low electron density, rie < lO'^cm"'^. We crudely 
estimate that the AGN can provide sufficient UV photons 
to account for the [OIV] luminosity 400 pc from the nu- 
cleus, as long as UV photons are not absorbed by interven- 
ing material as they travel from the AGN to the dust shell. 
A much more detailed photo-ionization study is required 
to understand the excitation of the [OIV] and [NeV] emis- 
sion. Previous reports of an ionization cone in Cen A were 
based on near-infrared ima ging of the central few arcseconds 
(|Brvant fc HunsteadllT999l '). Subsequent studies interpreted 
line emission in terms of a disk rather than ionization cone 
(Schrcicr ct al. 1998; Krajnovic ct al. 2005). Unfortunately 
the morphology of the near-infrared images is strongly af- 
fected by extinction and the warp disk models have not been 
good enough to accurately predict the extinction in the cen- 
tral few arcsecond of the nucleus. Cen A might be the only 
active galaxy in which mid-infrared spectroscopy has found 
evidence for high ionization species such as NeV at hundred 
pc distances from the nucleus. As far as we know Cen A 
hosts the only ionization cone that has been resolved with 
observations from the Spitzer Space Telescope. 

We see evidence for warm molecular hydrogen coinci- 
dent with the peak in [OIV] in the odd pure rotational odd 
transitions S(3) and S(5). The S(7) and S(2) transitions are 
also detected but at weaker levels. The S(3) and S(5) emis- 
sion also lies in the vicinity of the dust shell that is most 
prominent in the 11.3/im dust (PAH) emission feature. A 
two temperature component model can fit the rotational 
line ratios and implies that there is warm molecular hydro- 
gen with temperatures in the range 250-720K. The temper- 
atures are warmer than seen in nucleus of non-active nearby 
galaxies, similar to those of LINERS and Seyferts but lower 
than exhibited by supernova remnants ULIRGS (as com- 
pared to studies bv lRoussel et al.ll2007l , Iffigdon etai1l2006l 



8 



and lNeufeld et al.ll2007l ) . Near the jet axis, the column depth 
of warm molecular hydrogen is N{H2) ~ 10^°cm-^ similar 
to that estimated from the infrared co ntinuum emission of 
the dust shell by Qui llen et alj (|2006al ). This suggests that 
gas associated with the dust shell has been heated near the 
jet axis. There is probably more than one excitation process 
as [OIV] and [NeV] emission are detected primarily along 
the jet axis whereas the higher S pure rotational molecular 
hydrogen lines are detected there and in the vicinity of the 
dust shell. 

Previous studies of the pure-rotatio nal molecular hydro- 
gen lines in extra-ga l actic objects ( e.g., iHigdon et all 1200' 
iPanuzzo et al] l2007l : iRoussel et alll2007l : lOgle et all m& 
have not well resolved the emission. The association of the 
warm molecular hydrogen gas with a shell is most similar to 
phenomena exhibited by Galactic supernova remna nts where 
the b last wave encounters molecular clouds (Neu feld et al.l 
I2OO7I ). The physical conditions estimated from the molec- 
ular hydrogen observations are similar in properties (col- 
umn depth, temperatures and fraction of gas in the two 
temperature components) to the parameters estimated by 
[Neufeld ct al. (2007) in 4 Galactic supernova remnants. This 
suggests that theory of interstellar sho ck waves be apphed to 
interpreting observations of the shell. iNeufeld et al.l (|2007l ) 
associates the warmer molecular hydrogen component re- 
sponsible for the higher S pure rotational transitions with 
disso ciative shocks. These re quire shock velocities > 70km/s 
fe.g.. lDraine fc McKedl 19931 ) . A physical scenario and model 
accounting for the shell's structure and energetics is cur- 
rently lacking. Deep optical spectroscopic and radio studies 
are particularly needed to better constrain gas motions and 
physical conditions in this shell. 
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Figure 1. Continuum emission from the LH spectral cube from 
OAfim wide bands. Contours are evenly spaced with the lowest 
contour and spacing at 0.01 MJy/SR. a) Continuum centered at 
25.0/xm. b) Continuum centered at 32.0/im. 




Right Ascension 

Figure 2. Dust emission feature at 11.3/im. No continuum has 
been subtracted as the emission feature dominates the spectrum 
by a factor of 3-8. The minimum contour and spacing is approx- 
imately 10 MJy/SR in the peak of the line. The black contours 
show t he oval dust shell previously described by lOuillen et al.l 
l|2006al ). 
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Figure 3. Continuum subtracted line emission images in the LH 
spectral cube, showing emission in the folded star forming disk. 
Contours are evenly spaced. The lowest contours and spacing for 
the [SlII](33.48lAini) and [Sill] (34.815m) images are 1.0 and 1.5 
xlO"*crg cm-2 g-i SR-^, respectively, a) For [SIII](33.481/Ltm). 
b) For [SiII](34.815/im). 



Figure 4. Continuum subtracted line emission images in the LH 
spectral cube, showing emission in the folded star forming disk. 

Contours are evenly sijaccd. The lowest contours and spacings for 
[Fcll](25.988^m) and [Felll] (22.925/ini) images with lowest con- 
tour arc 0.05 and 0.025 xlO^**erg cm^'^ s^^ SR^^, respectively, 
a) For [FeII](25.988/im). b) for [FeIII](22.925/im). 



Figure 5. Continuum subtracted line emission images in the LH 
spectral cube, showing emission in the folded star forming disk, 
for the H2S(0)(28.22/im) line with lowest contour at and spacing 
at 0.01 xlO-^erg cm-^ s'^ SR-^. 
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Figure 7. Radio emission at 5Ghz shown as grayscale with 
[OIV] {25.9/^m) contours . The radio map i s a 6" resolution map of 
the inner radio lobes by iHardcastle et alj |20oi). The [OIV]25.9/i 
and [NeV]24.3/im line emission are oriented approximately but 
not exactly along the jet axis. 



Right Aicanstqn 

Figure 6. Continuum subtracted line emission images in the 
LH spectral cube, showing emission near the jet axis. Contours 
are evenly spaced. The lowest contours and spacing for the 
[OIV]{25.890Atm) and [NeV]{24.318Atm) images are 0.5 and 0.1 
xlO-^erg cm-2 s"! SR-\ respectively, a) For [OIV] (25.890Atm). 
b) For [NeV](24.318/im). 
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Figure 8. Continuum subtracted line emission images from 
the SL spectral cube, showing the H2S(3)(9.665/im) and 
H2S(5)(6.909^im) lines with lowest contours at 1 xlO~^erg cm~^ 
SR~^. The higher S rotational lines from H2 exhibit different 
morphology than the H2S(0)28/im line that was seen primarily 
in the folded star forming disk. Emission in the higher S lines is 
seen above the disk. Contours are evenly spaced. Lowest contour 
and spacing is X 10"'' erg cm~-^ SR~^ for both lines. For a) 
H2S(3)(9.665Mm). For b) H2S(5)(6.909Mm). 



Figure 9. a) Emission in H2S(3){9.665/im) shown as gray scale 
ovcrlayed with [OIV] (25.890^) contours. Contours are evenly 
spaced with lowest contours at are 0.5 xlO~*erg cm~-^ s~^ SR~^ 
and a spacing 4 times this. The gray scale range for the H2S(3) 



image is O(white) to 7(black) X 10 "erg < 



SR-l. b) Same 



as a) except the H2S(3) emission is overlayed with contours of the 
ll.S/xm PAH dust emission feature (as shown in figure[2]l- Contour 
spacing and lowest level is 20 MJy/SR. 
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Figure 10. Flux at 7.7 and 8.6/im showing dust emission fea- 
tures. Contours are evenly spaced. No continuum subtra<;tion has 
been done. The star forming disk is evident as the parallelogram 
shaped feature corresponding to the folded disk. The dust shell 

is seen above and below the parallelogram feature. The minimum 
contour and spacing is approximately 10 MJy/SR in the peak of 
the line, a) The 7.7/im dust emission feature, b) The 8.6;im dust 
emission feature. 




Figure 11. Line emission in [NeII](12.81(Um). No continuum has 
been subtracted as the [Nell] line dominates the continuum by 
a factor greater than 2 everywhere. Contours are evenly spaced 
with lowest contour at 4 xlO~^erg cm~^ SR~^. 
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Figure 12. a) The dust emission feature at 7.7ixm. divided by 
that at ll.S^im. The lowest and highest contours are shown at a 
ratios of 0.5 (in shell) and 1.0 (in pajrallelogram) , with spacing 
of 0.1. Black is 1.0, white is 0.5. The 7.7 to 11.0//m dust feature 

ratio varies with the dust shell having the lowest ratio, b) The 
[Nell] (12.8/im) line divided by the 11/im dust emission feature. 
The lowest and highest contours are shown at ratios of 1.9 (in 
shell) and 3.0 (in parallelogram), with spacing of 0.1. Black is 
3.0, white is 1.9. The strength of the Nell] (12.8/im) compared to 
the 11.3/im dust feature also varies with the dust shell having the 
lowest ratios. 
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Figure 13. a) Spectrum from the star forming parallelogram. 
D refers to a dust emission feature, b) Spectrum from the jet re- 
gion. The pure rotational molecular hydrogen lines S(2)-S(7) are 
labeled as S2-S7. Molecular hydrogen lines are more prominent in 
the jet region than in the star forming parallelogram. In the star 
forming disk the [Aril] (6.985/im) line is brighter than the nearby 
pure rotational molecular hydrogen 8(2) line at 6.909/im whereas 
the opposite is true in the jet region. There is a change in the 
ratio of the dust emission features in these two spectra and in the 
ratio of the [NeII]/11.3 dust emission feature. 
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Table 1. Lines Detected 



Line 


Rest Wavelength(/im) 


[Fell] 


5.340 


H2S(7)J=9-7 


5.511 


Dust 


5.7 


Dust 


6.2 


H2S(5)J=7-5 


6.909 


r A ^TTl 

[AriiJ 


6.985 


Dust 


7.7 


Dust 


8.6 


XJ O /'Q^ T c; Q 
J^2'5(.'J )o=o-o 


9.665 


[SIV] 


10.51 


Dust 


1 1 Q 
ii.O 


Dust 


12.0 


H2S(2)J=4-2 


12.279 


Dust 


12.7 


[Nell] 


12.81 


Dust 


14.2 


[Felll] 


22.925 


[NeV] 


24.318 


[OIV] 


25.890 


[Fell] 


25.988 


H2S(0)J=2-0 


28.219 


[SIII] 


33.481 


[Sill] 


34.815 



The spectra in the star forming disk and jet region have the same shap e near 14.2;jm sug - 
gesting that the feature at that wavelength is a dust emission feature jSmith et al.ll200^ 
rather than due to emission from [NeV] (14.3/im). The H2S{4)J=6-4 line at 8.025/im is not 
detected. The emission at this wavelength is dominated by nearby bright 7.7 and 8.6/im 
dust emission features. The S(6) line at 6.1086/im is also not detected as it would have been 
overpowered by the dust emission feature at 6.2/xm. The even J molecular hydrogen rotation 
quantum states are singlets rather than triplets and so should be similar to 3 times fainter 
than nearby odd J transitions. At the [OIV] and H2S(3) peak south-west of the nucleus the 
[Aril] (6.985^m) is weaker than the H2S(5)(6.909/im) line. The opposite is true in the star 
forming disk. [FeII](5.3/im is detected in the disk and at the OIV peak near the jet axis. 
The line in the middle above list denotes the break in wavelength between the SL and LH 
wavelength coverage. 

Table 2. Nebular Line Fluxes at [OIV] peak 



Line 


Flux 


[OIV](25.9Mm) 


1.5 


[NeII](12.8Atm) 


9.0 


[NeV](24.3Mm) 


0.15 


[NeV](14.3AJm) 


< 0.5 


[SIII](33.5Atm) 


1.6 


[SiII](34.8/im) 


4.0 


[SIV](10.5Mm) 


1.6 


[FeII](26.0Mm) 


0.2 


[FeIII](22.9Atm) 


0.05 


[ArII](6.98Mm) 


~0.1 



The fluxes were measured in a square region centered at RA=13:25:25.7 DEC=-43:01:24 and 
ll"wide. that is a peak in [OIV] and H2S(3) emission (see figure O. The spectra from this 
region are shown in Figures 1 14b and ll3b . Fluxes are given in units of 10~^erg cm~^s~^SR~^. 
The [NeV]{14.3Aim) flux is an upper limit. We did not detect [NcVI] at 7.642^m. 



Table 3. Pure rotational H2 Line Fluxes at [OIV] peak 



Line 


Flux 


S(0){28.2/xm) 


0.21 


S(2)(12.2/xm) 


1.1 


S(3)(9.7Aim) 


5.2 


S(4)(8.0AJm) 


<3 


S(5)(6.9Mm) 


2.7 


S(6)(6.lAjm) 


<6 


S(7)(5.5Mm) 


0.5 



Fluxes were measured in the same region as listed in Table [2] and in the spectra shown 
in Figures 114b and 113b . Fluxes are given in units of 10~^erg cm~^s~^SR~^. For the S(4) 
and S(6) lines we estimate upper limits. The upper limits are high because of dust emission 
features. 




13*^ ^rtrt ^^s 




13*^ ^rtrt ^^s 



